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2^ : Abstract 

^ ■ A method for solving the shell-model eigenproblem in a severely truncated 

space, spanned by properly selected correlated states obtained by partition- 
'M" . ing the full configuration space, is proposed. The method describes in a 

practically exact way the low energy spectroscopic properties of nuclei, as 
exemplified in schematic models. The applicability of the method to heavy 
nuclei as well as in contexts different from the nuclear shell model is stressed. 
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Complete shell-model calculations in the region of medium- and heavy-mass nuclei be- 
come rapidly prohibitive because of the extremely large configuration space dimensions. 
Although the increasing power of the present day computers, togheter with a clever us- 
age of the Lanczos algorithm |l|, makes feasible exact calculations in configuration spaces 
of impressive dimensions, the largest complete calculation in a full major shell have been 
done, as to our knowledge, in the mass region around ^^Ni. The extreme redundancy of 
the computed quantities is another problem: in the mare magnum of the output of these 
calculations, one is only interested in the comparison of a limited number of eigenvalues 
(usually the lowest ones) with the observed energy levels and in the identification of the 
most significant components of the corresponding wave functions. On the other hand, the 
truncation of the configuration space may not be an easy task, since taking into account the 
effect of the excluded space requires a corresponding renormalization of the interaction. 

A main road to physically significant reductions of the configuration space is the use 
of some kind of correlated basis; typical examples are the multi-step shell model 0, the 
broken-pair approximation or the chain-calculation method |^. The problem common 
to these methods which make use of correlated bases is the overcompleteness of the set of 
basis vectors used in the calculation. Since this redundancy gives rise to spurious admix- 
tures, it is necessary to resort to special techniques to get rid of the spurious components, 
e.g. computing and analyzing the overlap matrix of the basis, a task which is lenghty and 
involved. 

In this letter we propose a method which accounts in a practically exact way for all the 
configurations needed to the description of the lower energy states, and is free of all the 
illnesses of the mentioned truncation schemes. The idea consists in describing the system of 
interacting nucleons in terms of correlated subsystems defined in orthogonal subspaces. Let 
us consider a system of A^ valence particles in a given model space, defined by a projection 
operator 

I = Y.\aiN){aiNl 
I 

where the ai refer to all possible independent ways the N particles can be distributed over 
the single-particle (s.p.) levels of the model space. We now partition the s.p. states included 
in the model space in two groups whose corresponding configuration spaces we call Pi(A^i) 
and P2{N2), so that 

^ = E E I P^'^N,) I pfN^) {pfN2 I iP^'^N, I, 

N1N2 ij 

where the quantum numbers Pi classify all possible ways of distributing Ai and A'2 particles, 
with A^i + A^2 = A^, on the s.p. levels of the two partitions, respectively. We accordingly 
separate the Hamiltonian of the system as 



H = Hi + H2 + H 



12, 



where Hi acts only on the Pi spaces, H2 on the P2 spaces and H12 is the interaction term 
between Pi and P2. Solving the eigenvalue problems for all the allowed values of A^i and A2 
within each partition separately. 

Hi I Ai A^i) = E^,{Ni) I Ai Ai), H2 \ X2 N2) = E^,{N2) \ X2 A2), 



allows to write 
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Eq.(|l|) is our main point. Since the wave function of the A^-particle system is written as a 
direct product of eigenfunctions of correlated subsystems defined in orthogonal spaces, none 
of the previously mentioned redundancy problems arises, i.e. it combines the advantages 
of using an orthonormal basis with those offered by a description in terms of correlated 
subsystems. Shell-model basis is not well suited to further reductions of its dimensions 
essentially because the residual interaction can be so strong that even configurations which 
are quite distant in energy from each other may be equally important in the description of 
a physical state. This is not the case when using correlated bases, where a significant part 
of the effect of the interaction is already included and even drastic truncations of the basis 
can be meaningful. As a consequence, an energy criterion to truncate expansion ([^) works 
well: only those basis states which are not too different in energy from the physical state 
one wants to describe are relevant. It is worth noting that, as is explicitely shown in the 
following example, there are no particular difficulties in computing the matrix elements of 
H12 in the basis (|l|), since they can be expressed in terms of quantities computed in each 
partition separately. 

As a concrete example, let us consider a general two-body shell-model Hamiltonian for 
a system of identical particles: 
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where Ni = Em aL«im . and A^^^j^^hk) = Emim2(^i"^i^2^2 I ^o^o) al^^al^^. 
The partitions yield the two Hamiltonians 



Hi = ^e^iVi + - J2 Gj^{iii2i3i4) A^j„Mo(^i^2) A 
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H2 = Y1 ^J^'J + 4 51 G'jo(iii2j3J4) ^joMoO'iia) Aj^Moihu) 
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and the interaction term 

^12 = T J2 '^Jo(^1^2ilj2) ^JoMo (^1^2) ^JoA/o (ilia) 
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where ai^ = (— l)'^™ai_m and Jq = (2Jo + 1)^/^. The F matrix elements are defined as 

Fj^ilahUd) = - E (2-/1 + 1) Gj^{Udch) I J" J^ ;^° I , (6) 

and Bj^j^,j^{lil2) are the particle-hole operators: 

BjoMoiW = E {hmil2m2 \ JqMo) al^^di^^^. (7) 

m\m2 

Let I A^i7i JiMi) be the eigenstates of Hi for a system of A^i particles distributed over the 
s.p. levels of the partition 1 of the model space, with corresponding eigenvalues E^-^j-^{Ni), 
and I N2J2J2M2) the eigenstates of H2 for a system of N2 particles distributed over the s.p. 
levels of the partition 2, with corresponding eigenvalues E^^j^{N2). 

A complete basis for the configuration space of A^ identical particles, distributed over 
all the s.p. levels of the model space, with total angular momentum JM, can be therefore 
written as 

I iVi, 71, Ji,iV2, 72,^2, iVJM) = E {J1M1J2M2 I JM) I NaiJiMi) I N2I2J2M2). (8) 

M1M2 

where A^i and N2 run over all the possible values such that A'^i + A^2 = ^• 

To exemplify the quality of the results obtainable when an energy truncation criterion 
is used, we consider for simplicity a pairing Hamiltonian 

H = ^e,Ni + W Goil, n 4(1) AS'), (9) 

I ^ w 

where we use the shorthand Gq{1-, V) = Gq{111'1') and Al{l) = Aqq{11), and restrict ourselves 
to a seniority zero approximation for the configuration space. The interaction term @ 
becomes 

H12 = lT.Goihj) Al{^) A,{j) + h.c.. (10) 

The basis for an A^-particle system is now | A^i7i, A^272; N) = \ Niji) \ N2'~f2)- The matrix 
elements of H between those basis states are: 



{N[yi,N^i,;N\ H I iVi7i,iV272;iV) = [E,,{Ni) + E,,{N2)\6^>N^6^>N^Sy^,^Sy^^^ (11) 



+-j:gs3) 



+ -^i7j7j(^l) -^^^7^72(^2) ^N^{N[+2) ^N!^{N^+2) 

where Xj^^ yiNi) = {Ni 7^ || Aq(/c) || {Ni — 2) 7;') are two-particle transfer amplitudes. 

As already pointed out, matrix elements of H12 are quite simple in structure and can be 

evaluated using quantities defined in each subspace separately. It is worth noting that this 



is not peculiar of a pairing interaction: as can be easily seen from the structure of Eq.(^), 
matrix elements of if 12 in the basis (§) can be written in terms of the the two-particles 
transfer amplitudes Xi^i^j^^^j^y^j'^{Ni) = {Ni,ji,Ji \\ A^j^ihk) \\ Nl,'y'i,J'i), of the particle- 
hole matrix elements {Ni^'^i^Ji \\ Bj^{lil2) \\ Ni,'y'i,J[) and of the one-particle transfer 
amplitudes (A^/,7i, Ji || a|^ || A^/,7;, J/), computed in each partition separately. 

To reduce the dimension of the eigenvalue problem we retain only the basis vectors 
I A''i7i, A^272; N) corresponding to values of [E^-^{Ni) + E^^{N2)] up to a truncation value E*: 
progressively increasing E*, one can see how the energies in the full space approximate the 
exact ones. We have chosen a model space of ten s.p. levels divided in two partitions of five 
s.p. levels. In Table I are reported the first three energies obtained with various truncations 
of the basis with A^ = 30 nucleons. The ten s.p. levels are O/7/2, O/5/2, 1^3/2, lpi/21 05^9/2, 
which define the partition 1, and Og'7/2, 1^5/2, l'^3/2, 2si/2, 0/iii/2, which define the partition 
2. The corresponding s.p. energies are, in MeV, 0.0, 3.0, 3.5, 4.0, 5.0, 13.0, 13.5, 15.0, 
15.5, 16.0 and we use a constant pairing strength G = 0.16 MeV, where Go{ll') = —ll'G. 
We obtain practically exact results already diagonalizing matrices of order 150 while the 
dimension of the full space basis is 59, 702. 

Just to show that the existence of a gap in the s.p. energies is not essential to the 
quality of the approximation, we have also considered the simple, although not at all trivial, 
example of a pairing problem in a model space of 32 equispaced doubly-degenerate levels. 
Exact solutions of this problem in the case of 32 particles have been found in Ref. for the 
energies of the first three states. We arbitrarily choose as partition 1 the first 16 levels and as 
partition 2 the remaining ones. We have generated almost exact solutions for each value of A^i 
and N2 within each subspace, with the corresponding two-particle transfer amplitudes, using 
the method of Refs. Q |^. The results, reported in Table II, show an impressive convergence 
towards the exact values: the best computed energies differ by less than 0.1% from the exact 
ones, while the corresponding excitation energies, 3.060 and 4.900, are practically exact. It 
is worth noting that we use a number of states which is an exceedingly small fraction of the 
full basis dimension 6.01 x 10*. These results are significantly better than those of Ref. [|^, 
where a careful analysis of truncations of the configuration space is made; this fact not only 
confirms the effectiveness of the present method but also suggests that it can be a viable 
technique to obtain almost exact solutions for other physically interesting systems like those 
encountered in the modelization of ultra-small superconducting grains ^ [^ . 

We have considered in the preceding discussion only two partitions, just to be definite. 
One can have, however, situations where a basis generated by multiple partitioning the 
configuration space would be required. Althoug we do not treat these cases explicitely in 
this paper, it is worth noting that the generalization of the method presented here to systems 
where multiple partitioning is needed is quite straightforward. The method can be naturally 
applied to systems of neutrons and protons, which can be conveniently treated working out 
the eigenproblems for neutrons and protons separately by further partitioning of their model 
spaces, and then diagonahzing the n-p interaction in the product basis; it also constitutes 
a natural framework to treat excitations outside a major shell, e.g. in the study of ground 
state correlations. 
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TABLES 

TABLE L Absolute energies, in MeV, of the first three v = states, N = 30, for the model 
space of ten s.p. levels. The first row gives the number of states to which the full basis is truncated. 

25 60 150 230 430 810 

'Yq 90.256 90.244 90.241 90.241 90.241 90.241 

El 104.286 104.065 103.990 103.972 103.962 103.960 

E2 106.474 106.058 106.010 105.996 105.991 105.989 



TABLE IL Absolute energies, in MeV, of the first three states for the model space of 32 

twofold-degenerate equispaced levels, A^ = 32 particles, level spacing 1 MeV, pairing constant 
G = 0.345 MeV. The first row gives the number of states to which the full basis is truncated. 

1057 2078 3120 3454 3999 exact 

'Yo 263.587 263.453 263.390 263.375 263.356 263.171 

El 266.599 266.498 266.443 266.434 266.416 266.278 

E2 268.496 268.352 268.292 268.271 268.256 268.071 



